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Abstract 

We present a statistical study of prominence and filament eruptions observed by the 
Atmospheric Imaging Assembly (AIA) aboard the Solar Dynamics Observatory (SDO). 
Several properties are recorded for 904 events that were culled from the Heliophysics Event 
Knowledgebase (HEK) and incorporated into an online catalog for general use. These 
characteristics include the hlament and eruption type, eruption symmetry and direction, 
apparent twisting and writhing motions, and the presence of vertical threads and coronal 
cavities. Associated flares and white-light coronal mass ejections (CME) are also recorded. 
Total rates are given for each property along with how they differ among hlament types. 
We also examine the kinematics of 106 limb events to characterize the distinct slow- and 
fast-rise phases often exhibited by hlament eruptions. The average fast-rise onset height, 
slow-rise duration, slow-rise velocity, maximum held-of-view (FOV) velocity, and maximum 
FOV acceleration are 83 Mm, 4.4 hours, 2.1 km s"^, 106 km s"^, and 111 m s"^, respectively. 
All parameters exhibit lognormal probability distributions similar to that of CME speeds. A 
positive correlation between latitude and fast-rise onset height is found, which we attribute 
to a corresponding negative correlation in the average vertical magnetic held gradient, or 
decay index, estimated from potential held source surface (PFSS) extrapolations. We also 
hnd the decay index at the fast-rise onset point to be 1.1 on average, consistent with 
the critical instability threshold theorized for straight current channels. Finally, we explore 
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relationships between the derived kinematics properties and apparent twisting motions. We 
hnd that events with evident twist have signihcantly faster CME speeds and signihcantly 
lower fast-rise onset heights, snggesting relationships between these values and flux rope 
helicity. 

Keywords: Prominences, Dynamics; Coronal Mass Ejections, Low Corona Signatures; 
Corona, Structures 


1. Introduction 

Prominences are among the most common and well-studied features of the solar atmo¬ 
sphere, having been chronicled since the Middle Ages and identihed as “cloud formations” 
in the corona in the 1850s (Tandberg-Hanssen, 1998). They are now recognized to be thin 
channels of relatively cool, dense plasma suspended in the corona by highly sheared or 
twisted magnetic helds above polarity inversion lines in the photospheric magnetic flux 
distribution. Bright in Ha emission when protruding from the limb, prominences appear 
dark compared to the surrounding chromosphere when seen on the disk, where they are 
referred to as hlaments. We will use the terms prominence and filament interchangeably. 
Readers seeking detailed background information may direct their attention to the re¬ 
cent book. Solar Prominences (Vial and Engvold, 2015), from which we reference several 
chapters, and the excellent review of prominence observations by Parenti (2014). 

This paper focuses on prominence eruptions observed by the Atmospheric Imaging As¬ 
sembly (AIA, Lemen et al. 2012) aboard the Solar Dynamics Observatory (SDO, Pesnell, 
Thompson, and Chamberlin 2012), launched in 2010. The AIA provides continuous, full- 
sun coverage in 7 extreme ultraviolet (EUV) and 3 UV channels with high time (12s) and 
spatial (0.6'') resolution. A number of case studies have been published on AIA observations 
of eruptive prominences. Sterling, Moore, and Freeland (2011), Tripathi et al. (2013), and 
Chen, Bastian, and Gary (2014) present separate observations of active region hlament 
eruptions, each found to be consistent with an erupting flux-rope model. Su and van 
Ballegooijen (2012, 2013) present observations and flux-rope modeling of a polar crown 
eruption, and Thompson (2013) further examines signs of twist in the same event. Williams, 
Baker, and van Driel-Gesztelyi (2013) estimate the mass of erupted hlament material in 
a particularly stunning event, and several authors have combined the SDO with the Solar 
Terrestrial Relations Observatory (STEREO, Howard et al. 2008) spacecraft to exa mi ne 
prominence eruptions from three separate viewing angles (e.g. Chifu et al. 2012; Koleva 
et al. 2012). Su, Lu, and van Ballegooijen (2012) present a survey of 45 quiescent promi¬ 
nence eruptions observed by the SDO that serves as a pilot study for our work, with results 
generally consistent with those to be described here. 

We have recorded some basic observational properties for 904 hlament eruptions observed 
by the AIA, compiled an online catalog designed to aid the community in identifying 
promising events for future research, and performed a kinematics study using events se¬ 
lected from the catalog. §2 describes the catalog and the observations it provides. §3 details 
the recorded properties, including some background information on their signihcance and 
a report of our results. These include the hlament type (§3.1), eruption symmetry (§3.2) 
and direction (§3.3), apparent twisting (§3.4) and writhing (§3.5) motions, the presence of 
vertical threads (§3.6) and coronal cavities (§3.7), and white-light CME (§3.8) associations. 
§4 details our kinematics study. This involves height-time measurements of 106 limb erup¬ 
tions, which are ht with an analytic approximation to provide statistics on the slow- and 
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Figure 1. Start time (left) and position (right) data for filament eruptions found in the HEK and 
incorporated in the catalog. 


fast-rise phases often exhibited by hlament eruptions. Our methods are detailed in §4.1, 
the analytic approximation is described in §4.2, the impact of solar rotation is considered 
in §4.3, and our results are described in §4.4. §5 discusses interesting patterns that emerge, 
including a correlation between latitude and the fast-rise onset height (§5.1) along with 
suggestive relationships between apparent twisting motions and the derived kinematics 
values (§5.2). Finally, §6 summarizes our work. 


2. Online Catalog 

The catalog includes 904 hlament eruptions observed by the SDO between June of 2010 
and September of 2014. Times and positions are obtained from the Heliophysics Event 
Knowledgebase^ (HEK, Hurlburt et al. 2012), which incorporates a wide variety of phe¬ 
nomena and relies on reports from individual observers for hlament eruptions. Figure 1 
shows the event time and position distributions. An enhancement of eruptions near the 
limb is apparent, which is primarily due to the contribution of events behind the limb. We 
also note the general position of each eruption, hnding 401 disk events, 363 limb events, 
and 139 events originating behind the limb. 

For each event, we provide cutout and full-disk movies in the 171, 193, and 304 A AIA 
passbands. If available, full-disk observations are also provided from the 304 A channels of 
the Extreme Ultraviolet Imagers (EUVI, Wuelser et al. 2004) aboard the STEREO space¬ 
craft. Images are obtained primarily from the Helioviewer^ application program interface 
(API), which provides full-resolution JPEG2000 images that we then process separately to 
produce our content. In most cases, the 171 A images are processed from the original level 
1 data hies using a radial hlter described by Masson et al. (2014). All movies are provided 
for download and may also be played via YouTube. 

Users may sort and hlter the catalog based on the characteristics described in §3. Ad¬ 
ditional comments on each event from the HEK and catalog observers are also provided 
and may be used to search the catalog for specihc terms. Further details on the interface 


^Heliophysics Event Knowledgebase: http://www.lmsal.com/hek/ 
^Helioviewer: http://helioviewer.org/ 
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can be found at the catalog website^. Each event is also assigned an arbitrary rating based 
on how interesting a particular observer found it to be, which is intended to aid users 
who wish to browse through the large number of eruptions. A collection of movies from 
the highest-rated events is available as a YouTube playlist or as a downloadable TAR hie, 
which may also be useful for education and public outreach. The catalog table may also 
be downloaded as an IDL save hie or as a tab-delimited text hie. Finally, we note that the 
structure of the website is essentially the same as that of the Hinode/XRT Flare Catalog'^ 
and the Hinode & SDO Sigmoid Catalog® (Savcheva et a/., 2014). 

3. Statistical Properties 

The following subsections detail the properties recorded for each event, including brief 
background information on their signihcance and a report of our results. Each determina¬ 
tion was made through manual inspection by one of three observers with the exception of 
hare associations, which are automatically populated using data from the “Latest Events” 
archive in SolarSoft® and are not discussed further. Table 1 lists the event counts for each 
category, while Figures 2-9 provide representative examples and bar plots of the category 
distributions. All images are single frames from the event movies provided in the catalog 
with no additional processing. Clicking the ID numbers listed in the electronic text (e.g. 
No. 0962) will launch the catalog website, where the associated movies can be found. 

3.1. Type 

We assign types based on the commonly-used, location-based scheme. Active region (AR) 
hlaments are found within active regions, intermediate prominences (IP) are found adjacent 
to or between active regions, quiescent (QS) hlaments are found well away from active 
regions in the quiet sun, and polar crown (PC) prominences are found at high latitudes. 
Figure 2 shows representative examples of each type. Note that PC may refer specihcally 
to hlaments that form at the boundaries of the polar coronal holes, but we use the term 
more generally to refer to any hlament with a latitude greater than ~50° N/S. Including 
ambiguous events (e.g. “AR?”), our sample includes 172 AR, 124 IP, 451 QS, and 132 PC 
hlament eruptions. An additional 20 eruptions are of indeterminable type. These events 
are generally behind the limb without available STEREO observations, though a few are 
very compact hlaments in small emerging hux regions that do not ht neatly into the four 
main groups. 5 more are classihed as transequatorial (TE) hlaments, the longest and rarest 
variety. Neither of these latter two groups are included in the discussions below, but we hnd 
that the TE events have characteristics most similar to QS eruptions. Other classihcation 
schemes exist based on morphology, dynamics, spectroscopic characteristics, and magnetic 
conhguration. A review of these, along with historical classihcations, is provided by Engvold 
(2015). It would be worthwhile to investigate a subset of the events and characteristics 
described below using additional type sets, particularly the scheme based on photospheric 
hux distribution advanced by Mackay, Gaizauskas, and Yeates (2008). 

The extent to which each hlament erupts is also categorized, guided by conventions 
described by Gilbert, Alexander, and Liu (2007). A full (F) eruption is one in which 


^SDO Filament Eruption Catalog: http://aia.cfa.harvard.edu/filament/ 
'^Hinode/XRT Flare Catalog: http://xrt.cfa.harvard.edu/flare_catalog/ 
^Hinode & SDO Sigmoid Catalog: http://aia.cfa.harvard.edu/sigmoid/ 
®SolarSoft Latest Events: http://www.lmsal.com/solarsoft/latest_events/ 


SOLA: filament_catalog_paper.tex; 11 May 2015; 0:17; p. 4 




Prominence Eruptions Observed by the SDO 


Table 1. Properties and event counts. Ambiguous (e.g. “Y?”) event counts are listed in parens. 

Type 


Category 

Bin 

All 

AR 

IP 

QS 

PC 

TE 

? 

Type 

904 

135 (37) 

76 (48) 

392 (59) 

122 (10) 

4 (1) 

20 


Disk 

401 

92 

68 

214 

19 

5 

3 

Position 

Limb 

363 

50 

48 

178 

82 

0 

5 


Behind 

139 

30 

8 

58 

31 

0 

12 


Full 

285 (73) 

56 (13) 

39 (7) 

143 (42) 

36 (9) 

3 

8 (2) 

Eruption 

Partial 

266 (99) 

38 (23) 

38 (17) 

131 (49) 

53 (7) 

2 

4 (3) 

Confined 

122 (56) 

33 (9) 

15 (7) 

57 (29) 

14 (11) 

0 

3 (0) 


Other“ 

18 (2) 

5 (1) 

1 

10 (1) 

2 

0 

0 


Sym 

307 (126) 

41 (30) 

49 (13) 

174 (62) 

34 (17) 

2 (1) 

7 (3) 

Symmetry 

Asym 

232 (114) 

49 (19) 

29 (16) 

116 (48) 

33 (27) 

1 (1) 

4 (3) 


? 

125 

33 

17 

51 

21 

0 

3 


Radial 

625 (51) 

95 (12) 

79 (3) 

338 (26) 

97 (7) 

2 (2) 

14 (1) 

Direction 

Non-R 

73 (35) 

11 (6) 

12 (4) 

39 (17) 

10 (6) 

0 

1 (2) 


Sideways 

72 (26) 

36 (8) 

18 (5) 

12 (8) 

5 (4) 

1 

0 (1) 


? 

22 

4 

3 

11 

3 

0 

1 


Yes 

408 (130) 

70 (31) 

65 (19) 

201 (63) 

58 (14) 

1 

13 (3) 

Twist 

No 

265 (68) 

44 (14) 

28 (8) 

143 (34) 

47 (10) 

1 (1) 

2 (1) 


? 

32 

13 

4 

9 

3 

2 

1 


Yes 

73 (71) 

26 (14) 

15 (10) 

27 (38) 

1 (7) 

0 

4 (2) 

Writhe 

No 

656 (61) 

102 (14) 

79 (13) 

350 (23) 

110 (10) 

5 

10 (1) 


? 

43 

16 

7 

13 

4 

0 

3 

Vertical 

Yes 

384 (58) 

10 (6) 

52 (16) 

246 (30) 

69 (6) 

4 

3 

Threads^ 

No 

296 (26) 

118 (8) 

42 (6) 

105 (11) 

26 

0 (1) 

5 


Yes 

97 (27) 

2 (2) 

10 

47 (12) 

37 (12) 

0 

1 (0) 

Cavity“ 

No 

165 (13) 

35 

25 

79 (8) 

25 (5) 

0 

1 


During 

47 (14) 

9 (2) 

9 (3) 

24 (8) 

3 

0 

2 (1) 

CME 

Yes 

609 (44) 

108 (12) 

85 (8) 

307 (18) 

91 (5) 

4 

14 (1) 

No 

225 (26) 

50 (2) 

26 (5) 

112 (14) 

32 (4) 

1 

4 (1) 

Flare 

Yes 

222 (15) 

94 (2) 

50 (3) 

68 (9) 

7 (1) 

1 

2 


No 

682 (15) 

78 (2) 

74 (3) 

383 (9) 

125 (1) 

4 

18 


“Indicates non-eruptive events that are excluded from all other counts. 
^Vertical threads counts exclude events behind the limb. 

“Cavity counts include limb events only. 
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Figure 2. Examples of the four main types: active region (AR, No. 0631), intermediate (IP, No. 0667), 
quiescent (QS, No. 0977), and polar crown (PC, No. 0621). 


the bulk of the fil am ent mass escapes the AIA £eld-of-view (FOV) along with the full 
magnetic structure, a partial (P) eruption is one in which the magnetic structure erupts 
either completely or partially, expelling some or none of the hlament mass, and a confined 
(C) or failed eruption is one in which none of the hlament mass nor magnetic structure 
escapes. It can be difficult in practice to apply these straightforward dehnitions because 
all events exhibit some draining, whereby expelled material hows back down the legs of 
an erupting prominence. This returned mass may constitute an apparently signihcant 
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Figure 3. Examples of symmetric (left, No. 0745) and asymmetric (middle, No. 0405) eruptions. The 
right panel shows the number of events found for each type, and shaded regions denote ambiguous events 

(e.g. “S?”). 


volumetric fraction of the material seen, but it is not necessarily clear what fraction of 
the actual mass escapes given that the volume of the structure is evolving and that the 
304 A observations we primarily rely on are optically thick even at low column densities. 
For this reason, we relax the >90% criterion used by Gilbert, Alexander, and Liu (2007) 
and dehne a full eruption as one in which the magnetic structure erupts completely and no 
more than a third of the hlament material appears to drain back down. Question marks 
are used to denote events near this boundary. 

It can also be difficult to make a determination for faint events that become too diffuse 
to be tracked to the edge of the FOV or distinguished from the background chromosphere 
on the disk; these are labeled as full if they appeared to be erupting fully prior to becoming 
untraceable. The very few instances in which a hlament cha nn el erupts with little or no 
hlament material contained within it are labeled as partial eruptions (e.g. No. 0478). 
Finally, note that these labels are of course based on the behavior of the hlament and 
associated magnetic structure to the extent to which they are visible in the three passbands 
we examine. It may be possible, for instance, for an apparently conhned event to produce a 
detectable CME either because an eruption was stimulated in the overlying held or because 
some part of the prominence structure escaped without producing a noticeable signature 
in our observations (see §3.8). Given these dehnitions and including ambiguous events (e.g. 
F?), the catalog includes 358 full, 365 partial, and 178 conhned eruptions. An additional 
20 events are labeled as other (O), which include non-eruptive dynamics such as hlament 
activation. These events are included in the online catalog but are excluded from all counts 
in this paper. 

3.2. Symmetry 

Filament eruptions are often divided into symmetric and asymmetric categories depending 
on whether they exhibit a loop-like geometry or one that favors the liftoh of a single foot- 
point. This distinction is particularly useful for studies of potential eruption mechanisms 
and magnetic conhgurations. For instance, Tripathi, Isobe, and Mason (2006) examined 
EUV brightenings associated with both types and found that in the symmetric case, 
brightenings propagated from the middle to each end, while in the asymmetric case, they 
moved from the erupting end toward the tethered one. Based on this, they suggest a three- 
dimensional extension to the standard two-dimensional flare model in which the EUV 
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Figure 4. Examples of non-radial (left, No. 0797) and sideways (middle, No. 0332) trajectories. Arrows 
denote the approximate trajectories. The right panel shows the total number of radial (R), sideways (S), 
and non-radial (NR) events, and the shaded regions denote ambiguous events (e.g. “NR?”). 

brightenings are caused by successive reconnection events. Liu, Alexander, and Gilbert 
(2009) further divided asymmetric events into “whipping” and “zipping” types for which 
the active leg either whips upward dramatically or moves along the polarity inversion line 
like the unfastening of a zipper, with each case exhibiting distinct hard X-ray patterns and 
implying different held conhgurations. 

We distinguish simply between the symmetric (S) and asymmetric (A) cases based 
primarily on inspection of 304 A observations. This distinction may be obscured by the 
orientation of the hlament and its two-dimensional projection, such as when the hlament 
has a predominantly east-west axis and is positioned on the limb. For this reason, there are 
a large number of ambiguous (e.g. “S?”) and indeterminable (“?”) events. Figure 3 provides 
examples of both cases along with a bar plot showing the frequency distribution. Including 
ambiguous events, we hnd that 48% of eruptions are symmetric, 38% are asymmetric, and 
14% are indeterminable. The respective percentages of symmetric and asymmetric events 
by hlament type are 41% and 40% for AR, 50% and 36% for IP, 52% and 36% for QS, and 
39% and 45% for PC eruptions. Symmetric eruptions are thus somewhat more common 
for each type except polar crown hlaments. It is possible that the longer length of PC hla- 
ments could provide more opportunities for one footpoint to be preferentially destabilized 
by neighboring events, however these diherences may not be very signihcant given large 
number of questionable and indeterminable events. More signihcantly, symmetric events 
are around 1.5 x more likely to be full eruptions than asymmetric ones. This is primarily 
because the asymmetric geometry allows mass to more easily drain out of the hlament 
before it is fully expelled and because asymmetric events frequently leave behind a section 
of the hlament at the end opposite to the initial liftoh. 

3.3. Direction 

The directions of hlament eruptions and associated CMEs are inhuenced by the surround¬ 
ing magnetic environment, which may dehect an eruption away from a simple radial 
trajectory. Such interactions between small- and large-scale magnetic structures are of 
particular interest in the space weather context because CME dehection inhuences an 
event’s potential to impact Earth (Kay, Opher, and Evans, 2014). Gopalswamy (2015) 
provides a recent review of non-radial prominence trajectories, which generally match that 
of the associated white-light CMEs (Simnett, 2000) but often exhibit a greater degree 
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of non-radiality (Panasenco et al, 2013). Eruptions will tend to travel along the paths of 
least resistance, deflecting away from regions of higher magnetic energy density (Gui et al, 
2011) such as coronal holes and active regions. 

We distinguish between three types of trajectories: radial (R), non-radial (NR), and 
sideways (S). Non-radial events are somewhat inclined from a radial trajectory and are 
generally deflected by large scale features such as coronal holes and pseudo-streamers. Side¬ 
ways events have a substantial tangential component initially and are generally deflected 
away from active regions, often adopting a radial trajectory while still within the AIA 
FOV. Figure 4 provides examples of the latter two cases along with a bar plot showing 
the category distributions. The non-radial example is deflected southward by an overlying 
pseudo-streamer and the sideways example is ejected tangentially out of its host AR. 
Including ambiguous events (e.g. “NR?”), we hnd that 75% of eruptions are radial, 12% 
are non-radial, 11% are sideways, and 2% are indeterminable. By type, the respective 
percentages of non-radial and sideways events are 10% and 26% for AR, 13% and 19% for 
IP, 12% and 4% for QS, and 12% and 7% for PC eruptions. Similar non-radial rates of 
around 12 ± 1% are found for full, partial, and conhned eruptions. However, these groups 
exhibit sideways trajectories at different rates of 7%, 12%, and 16%, respectively. Sideways 
eruptions are more frequently conhned or partial perhaps because they are more likely to 
encounter neighboring magnetic structures that might arrest the entire eruption or divert 
some ejecta along a different set of held lines. Alternatively, eruptions may be conhned due 
to the prevailing strength of the overlying held and recorded as sideways if the material is 
diverted tangentially before draining back to the surface. 

Future work on quantifying the extent to which eruptions deviate from the radial di¬ 
rection and on what prompts this deviation might start from the events identihed here. 
In particular, there is at present no clearly-articulated model describing the evolution of 
the magnetic held in the sideways ejection events from ARs. One hypothesis is that such 
trajectories arise from conhgurations in which the overlying coronal arcade in the center 
of the active region is too strong for the hux rope to break through radially. Instead, 
the system may hrst destabilize in the weaker-held region at the active region periphery. 
This might explain why hlament plasma is seen to move horizontally from strong to weak 
held regions at the start of some events. However, if the stronger-held side is suddenly 
diminished by hux cancellation or footpoint motions, the hlament may also eject toward 
the weakened stronger-held side. Thus, new insights might be gleaned from magnetic held 
modeling along with a careful examination of potential changes in the photospheric hux 
distribution for some events in our sideways sample. 

3.4. Twist 

Apparent twist in hlament eruptions is particularly intriguing because of its relationship 
to stored energy release and interactions with nearby magnetic structures. Several types 
of twist may be observed, each with diherent implications. For instance, the “roll ehect” 
refers to rotation about the hlament axis whereby the loop apex rolls to one side, creating 
twist of opposite signs in the two legs (Martin, 2003). Panasenco et al. (2011) explore 
this ehect for three eruptions observed by STEREO, and Panasenco et al. (2013) consider 
its connection to prominence trajectories, showing a relationship between roll and the 
presence of large-scale features like coronal holes and pseudo-streamers that may induce 
rolling motions through dehection. Murphy et al. (2012) suggest that rolling may also arise 
from an ohset between a rising hux rope and the corresponding CME current sheet. 
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Figure 5. Examples of events exhibiting twist (left: No. 0994, right: No. 0696). The right panel shows 
the number of events that do (Y) and do not (N) show twist, and the shaded regions denote ambiguous 
events (e.g. “Y”?). 


Su and van Ballegooijen (2013) model a PC ernption, reprodncing rolling motions ob¬ 
served at low heights throngh a continnons injection of twist from the active region side 
that creates an asymmetric held distribution, possibly inducing rolling motions through 
asymmetric reconnection similar to that described by Murphy et al. (2012). Su and van 
Ballegooijen also describe further rotating motions attributed to the relaxing of dipped 
held lines after the eruption (see also Thompson 2013). Romano, Contarino, and Zuccarello 
(2003) explore untwisting motions in detail, which can often be seen in one or both legs 
during eruptions (e.g. Yan et al. 2014). Mass how during an eruption may also reveal 
twisted structure inherent to the hlament. Finally, a distinction can be drawn between 
twist, which generally refers to the degree to which held lines are wound around a central 
axis, and writhe, the helical deformation (twist) of the axis itself that may lead to the 
kink instability (Torok, Berger, and Kliem, 2010). These properties may be combined into 
the mathematical measure helicity, which quantihes the degree to which a magnetic held 
is twisted, writhed, and linked. Helicity is positively related to free magnetic energy, and 
helicity accumulation is thought to be very important in leading to an eruption (Tziotziou 
et al, 2014). 

In this work, we assess simply whether (un)twisting motions of any type are apparent in 
each eruption, and references to “twisting motions” here and in later sections indicate the 
presence of at least one of the twist signatures described above. It was our initial intent to 
subdivide twist into categories related to the previous discussion, but we found that the 
careful inspection required to make these distinctions was untenable for our large sample 
and additional aims. Future work on this topic is planned, and some preliminary results 
were shown by McKillop et al. (2014). However, we do note that the untwisting motions 
in hlament legs and mass how that reveals twisted structure are the most commonly- 
observed of these phenomena in our events. Figure 5 shows two examples and the event 
count distribution. Including ambiguous events (e.g. “Y?”), we hnd that 60% of events 
show signs of twist, 37% do not, and the remainder are indeterminable. By type, the 
percentage of events with evident twist are 59% for AR, 68% for IP, 59% for QS, and 55% 
for PC eruptions. We therefore hnd that a majority of all hlament types exhibit signs of 
twist and that the likelihood is highest for IP eruptions. 

We also hnd that events with apparent twisting motions have faster CME speeds on 
average, and that IP events are generally the fastest of the four types (see §4.4 and 
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Figure 6. Examples of events exhibiting writhe (left: No. 0126, right: No. 0927). The right panel shows 
the number of events that do (Y) and do not (N) show writhing, and the shaded regions denote ambiguous 
events (e.g. “Y”?). 

§5.2). Full and partial eruptions exhibit twisting motions at similar rates of 61% and 
66%, respectively. Interestingly, confined eruptions are significantly less likely to exhibit 
signs of twist, with a rate of 44%, despite including a comparatively higher rate of events 
with writhing motions (see §3.5). We hypothesize that events exhibiting twisting motions 
will have had greater helicity on average, which implies a greater release of magnetic 
energy during the eruptions. The lower rate of twisting motions in co nfi ned events thus 
might be understood in terms of the energy available to push through the overlying held. 
Relationships between apparent twisting motions and eruption kinematics will be explored 
in §5.2. 


3.5. Writhe 

The helical kink instability is an ideal magnetohydro dynamic (MHD) instability that affects 
columns of plasma with strong axial currents, such as magnetic hux ropes, and is so 
named for the characteristic inverted gamma shape of kinked rope. As described in the 
previous section, the helical deformation of a hux rope’s axis that may ultimately produce 
a kinked shape is generally referred to as writhe (Torok, Berger, and Kliem, 2010). The 
kink instability is one of the principal mechanisms thought to drive prominence eruptions 
(see §5 and review by Fan 2015), and numerical simulations of kink-unstable hux ropes 
have been very successful in reproducing characteristics of observed events (Torok and 
Kliem, 2005), suggesting that coronal hux ropes exist prior to eruption (Gibson and Fan, 
2006). An overview of how the kink instability manifests in solar observations is given by 
Gilbert, Alexander, and Liu (2007), and a recent ehort to quantify the twist in a likely 
kink-unstable event is presented by Yan et al. (2014). 

We note the appearance of writhing based on inspection of the 304 and 171 A obser¬ 
vations. Figure 6 provides two examples along with a bar plot showing the distribution 
of events with and without. Including ambiguous events (e.g. “Y?”), we hnd that 79% of 
eruptions show no writhe, 16% do, and another 5% were indeterminable. However, fully half 
of the 144 writhed events are labeled as ambiguous. This is due to the difficulty in positively 
identifying the inverse-gamma structure, which may be obscured by the orientation or 
dynamics of the event. By type, we hnd that 23% of AR, 20% of IP, 14% of QS, and 6% of 
PG eruptions exhibit writhing motions. 62% and 23% of the writhed events are partial and 
conhned eruptions, respectively, which are signihcantly greater than the general population 
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Figure 7. Examples of vertical absorption threads on the limb (left, No. 0578) and disk (middle, No. 
0797). The right panel shows the number of events with (Y) and without (N) vertical threads, excluding 
events behind the limb, and the shaded regions denote ambiguous events (e.g. “Y?”). 

rates of 40% and 20%. This is consistent with expectations from Gilbert, Alexander, and 
Lin (2007), who note that kinking strnctnres are highly nnlikely to prodnce fnll ernptions 
becanse the hlament material needs to be concentrated at the center of the flnx rope and 
the ernption needs to be snfhciently fast to limit draining. We do hnd that 15% of onr 
kinked sample (22 events) are fnll ernptions, which might be greater than expected by 
Gilbert, Alexander, and Liu. This is largely due to our comparatively lax dehnition of 
“full” (see §3.1), which permits more draining, however at least a few events exhibiting 
writhing motions appear to be full eruptions even with the stricter dehnition (e.g. No. 0770 
& No. 0207). It is also important to note that a number of additional events classihed as 
having apparent twisting motions in the previous section might also have been inhuenced 
by the kink instability but had underdeveloped or obscured inverse-gamma structures. 

3.6. Vertical Threads 

Thin threads, revealed primarily by high-resolution Ha imaging, are the ubiquitous build¬ 
ing blocks of solar hlaments (Lin, Martin, and Engvold, 2008). These are held-aligned and 
oriented horizontally, parallel or somewhat inclined with respect to hlament spines and 
barbs. In addition to horizontal threads, vertical threads, perpendicular to the spine, are 
frequently observed in quiescent and polar crown “hedgerow” prominences (Engvold, 2015). 
Vertical threads are challenging to explain given the predominantly horizontal orientation 
of the local magnetic held (Paletou et al, 2001), and several scenarios have been proposed. 
For instance, van Ballegooijen and Granmer (2010) suggest that vertical threads form in 
a tangled held environment, Ghae (2010) suggests that they form from stacks of sagging 
horizontal threads, and Low et al. (2012a,b) suggest that vertical structures may arise 
from cross-held mass transport after a breakdown of the frozen-in magnetic held condition, 
perhaps in concert with magneto-thermal convection as described by Berger et al. (2011). 

We note the visibility of vertical threads (or more generally, threads perpendicular to 
the hlament spine) based primarily on inspection of the 193 A images. The threads appear 
dark in the EUV, as emission is absorbed and not produced by the neutral hydrogen and 
helium of the partially-ionized prominence plasma, which typically has HI column densities 
of ~10^® cm"^ (Labrosse et al, 2011). Figure 7 shows examples of vertical threads seen 
on the limb and on the disk, along with a bar plot showing the number of events with 
and without. Including the ambiguous events (e.g. “Y?”), and excluding those behind the 
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Figure 8. Examples of a coronal cavity (left, No. 0452) and a cavity that became apparent during the 
eruption (middle, see movie: No. 0953). The right panel shows the number of limb events with (Y), without 
(N), and during (D), and the shaded regions denote ambiguous events (e.g. “Y?”). 


limb, we find that 58% of events exhibit vertical threads and 42% do not. By type, the 
percentage of events with vertical threads are 11% for AR, 59% for IP, 70% for QS, and 74% 
for PC ernptions. The size of a filament infinences whether the threads themselves or their 
orientations will be distingnishable, and we find rates that grow with the characteristic size 
scales of the fonr filament types. Confined events are also mnch less likely to exhibit vertical 
threads (43%), perhaps dne to their generally smaller sizes. This may be cansed by limited 
spatial resolution, or it may be that compact filaments are indeed less likely to contain 
vertical threads. We also note that the 16 AR events that do exhibit vertical threads are 
generally more extended, with higher pre-emption heights, than average for that type. 
Several of these events are also on the boundary between what might be considered an 
“active region” versus “intermediate” filament or are found within decaying active regions 
with well-separated polarities. Thus, we find the presence of vertical threads to be the 
norm outside of active regions, and we note that our rates may be underestimated because 
some small or unfavorably-oriented filaments were likely excluded. 

3.7. Cavity 

Coronal cavities are elliptical depressions in the intensities above and surrounding promi¬ 
nences observed at optical, EUV, and soft X-ray wavelengths. They can be understood 
as the limb projections of filament channels on the disk, which may or may not contain 
a filament, and their interiors are less dense (Fuller and Gibson, 2009) and often hotter 
(Reeves et al., 2012) than the surrounding plasma. Cavities are commonly modeled as 
magnetic fiux ropes (e.g. Low and Hundhausen 1995; Fan 2012) and may be inherent 
components of filament systems that form from helicity accumulation (Mackay, 2015). At 
least a third of cavities erupt (Forland et al, 2013), though many of those that do not 
may simply erupt sometime after they were observed at the limb. CMEs often exhibit a 
three-part structure that includes a bright plasma pileup followed by a dark cavity and 
bright prominence core (Hundhausen, 1999). A recent review of cavities is provided by 
Gibson (2015). 

We note the visibility of associated cavities based on inspection of 193 A observations, 
which are sensitive to material around 1.2 MK in non-haring regions. Of the AIA passbands, 
cavities are most clearly delineated in this channel along with the 211 A band, which is not 
included in the catalog. No (N) indicates that a cavity was not apparent, yes (Y) indicates 
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Figure 9. Left: Distribution of events with and without associated white-light CMEs. Middle: Overall 
CME speed distribution. Right: Cumulative CME speed distributions by type. 


that a cavity was visible prior to eruption, and during (D) indicates that a cavity became 
visible during the eruption. Figure 8 provides two examples along with a bar plot of the 
results for all 363 limb events in the catalog. Including the ambiguous events (e.g. “Y?”), 
we hnd that 49% have no obvious cavity prior to eruption, 34% do, and 17% have cavities 
that become visible after the eruption onset. 

By hlament type, we hnd that cavities are apparent before and during the eruptions 
for 8% and 22% of AR, 23% and 25% of IP, 33% and 18% of QS, and 60% and 4% of 
PC eruptions, respectively. It is important to note that cavities are most apparent when 
the hlament is large, oriented along the line-of-sight, and free of neighboring structures, so 
higher rates are expected for PC and QS prominences. Conversely, the compactness and 
proximity to active region arcades of AR and IP prominences explains why associated cav¬ 
ities are less frequently seen prior to eruptions and more frequently seen during eruptions, 
when motion makes the coherent structures more apparent. Finally, we note that additional 
processing such as radial hltering (Forland et al., 2013) can improve cavity visibility, which 
may have resulted in additional detections if applied to our 193 A observations. 

3.8. Coronal Mass Ejections (CMEs) 

Coronal Mass Ejections (CMEs) are large expulsions of plasma and magnetic helds that 
escape into the heliosphere. A strong association between CMEs and hlament eruptions 
has been known for many years from high statistical associations (~70%; Munro et al. 
1979) and Ho emission in CME cores (Sheeley et al, 1980). As described in the previous 
section, CMEs often exhibit a three-part structure that includes the overlying prominence 
cavity and embedded hlament. A review of CME and prominence eruption associations is 
provided by Webb (2015). Complementary reviews on prominence eruption dynamics and 
space weather implications are also provided by Gopalswamy (2015) and Lugaz (2015), 
respectively. White-light coronagraphs have historically been the primary means by which 
CMEs are detected, and our associations are based on observations from the Large Angle 
Spectroscopic Coronagraph (LASCO; Brueckner et al. 1995) aboard the Solar and Helio¬ 
spheric Observatory (SOHO; Domingo, Fleck, and Poland 1995). Our dehnition of “CME” 
is thus instrument-dependent, but it is important to note that the term refers to a physical 
process that is not limited to coronagraph observations. 
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Table 2. CACTus CME Speeds 






CME Speed (km s"^) 

Category“ 

Events 

Min 

Max 

Mean ± SD^ 

Median ± MAD^ 

All CACTus 

6368 

91 

1950 

428 ± 277 

367 ± 182 

All Filament 

609 

41 

1840 

425 ± 228 

390 ± 158 

Active Region 

85 

140 

1320 

466 ± 219 

425 ± 163 

Intermediate 

52 

149 

1710 

515 ± 272 

449 ± 182 

Quiescent 

266 

41 

1840 

391 ± 204 

365 ± 143 

Polar Crown 

85 

107 

1220 

351 ± 184 

337 ± 133 

Full 

264 

112 

1740 

461 ± 225 

422 ± 157 

Partial 

210 

41 

1230 

390 ± 208 

361 ± 151 

Twist = Yes 

306 

108 

1740 

432 ± 233 

380 ± 160 

Twist = No 

156 

41 

1840 

380 ± 210 

374 ± 147 


“Subgroups exclude ambiguous events (e.g. Type = AR?). 

^SD = Standard Deviation & MAD = Median Absolute Deviation. 


We search for matching events using the Computer Aided CME Tracking catalog (CAC- 
Tus^), which automatically detects CMEs, measuring both velocity and angular size. The 
speeds are averages of linear hts to radial height-time prohles that span the angular width 
of the CME and extend from 1.5 to 30 R© in height above the limb (Robbrecht and 
Berghmans, 2004). These measurements reflect all parts of the CME, but note that there 
is variation within the three-part CME structure, with leading edges tending to be 30 to 
40% faster than core speeds (Gopalswamy et ai, 2003; Maricic, Vrsnak, and Rosa, 2009). 
In the few cases where LASCO or CACTus data are unavailable, we note the likelihood of 
a CME in the catalog with a question mark (“e.g. Y?”). We hnd that 72% of our eruptions 
have associated CMEs and 28% do not. Similar rates are found when the events are divided 
by type, with CME-associations found for 70% of AR, 75% of IP, 72% of QS, and 73% of 
PC eruptions. 

Basic statistics are listed in Table 2 for several groupings, and Figure 9 shows the CME 
speed distributions, which exhibit lognormal prohles that have been reported in other 
works (Yurchyshyn et al, 2005; Zhang and Dere, 2006; Bein et ai, 2011). We will see in 
§4 that this pattern emerges in the low corona. The average speed of CMEs associated 
with our hlament eruptions is 430 km s"^, which is essentially the same as for all CMEs 
detected by CACTus over the SDO mission lifetime. Moon et al. (2002) compare similarly- 
observed CMEs associated with hares versus hlament eruptions. They found a similar 
average velocity for hlament-eruption CMEs but also reported these to be ~9% slower 
than the general population, which is not found here. This discrepancy may be due to 
the inclusion of false-positives in the CACTus catalog from transients related to previous 
CMEs or intensity variations in the slow solar wind that bias the total CACTus distribution 
toward lower speeds (Robbrecht, Berghmans, and Van der Linden, 2009). 

We hnd, predictably, that full eruptions produce faster CMEs than partial eruptions by 
a factor of 1.2 on average. Full and partial eruptions have corresponding LASCO CMEs 
95% and 81% of the time, respectively. Those that do not produce CMEs are generally 


^CACTus CME catalog: http://sidc.oma.be/cactus/ 
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very small, perhaps becoming confined or too diffnse to be detected at larger heights, or 
are oriented on the disk such that the opening angle is unfavorable for viewing by the 
coronagraph. Conversely, confined events lack CMEs 96% of the time. The remaining 4% 
either stimulate an eruption in a higher or neighboring structure, or the escaping filament 
material and/or magnetic structure was not perceptible in the AIA observations. 

We also find that intermediate prominence eruptions are the most likely to be associated 
with CMEs and that their CMEs tend to be the fastest. On average, IP-CMEs are 1.1 x 
faster than AR-, 1.3x faster than QS-, and 1.5x faster than PC-CMEs. The cumulative 
probability distributions shown in the right panel of Figure 9 indicate that the AR and 
IP curves are quite close. To test if their difference is statistically significant, we apply the 
k-sample Anderson-Darling (A|) test, a non-parametric test that measures the deviation 
between an arbitrary number of empirical distribution functions (2 in this case) and pro¬ 
vides a significance level for which the null hypothesis that the groups are drawn from the 
same distribution can be rejected (Scholz and Stephens, 1987). The A| statistic for the IP 
and AR distributions is -0.50, corresponding to a p-value of 0.59, well above the 0.05 level 
below which the distributions can be considered significantly different. We conclude, then, 
that our finding of enhanced CME speeds for IP versus AR eruptions is suggestive but 
not statistically significant. We do find a significant relationship between apparent twisting 
motions and CME speed, which will be discussed in §5.2. 


4. Kinematics 

Filament eruptions often exhibit distinct slow- and fast-rise phases, which has been ex¬ 
plored in a number of case studies (e.g. Sterling and Moore 2004, 2005; Chifor et al. 2006; 
Sterling et al. 2007; Schrijver et al. 2008; Sterling, Moore, and Freeland 2011; Joshi and 
Srivastava 2011; Koleva et al. 2012). The slow portion is characterized by speeds of a few 
km s'^ or less that may persist for several hours or minutes for quiescent and active region 
filaments, respectively. Rapid acceleration to hundreds of km s"^ follows during the fast- 
rise phase, which may produce a CME or a confined eruption. The slow-rise phase is often 
attributed to gradual flux cancellation or emergence that builds and elevates the filament 
until it succumbs to an MHD instability or a fast-reconnection process that initiates the 
fast-rise phase (Sterling, Harra, and Moore, 2007; Sterling, Moore, and Freeland, 2011; Fan, 
2015). We will consider the underlying mechanisms in greater detail with respect to our 
results in §5. While previous low-corona studies have primarily relied on close inspection 
of individual events or a small ensemble, we examine the kinematics of 106 eruptions in 
our catalog to determine general properties of the slow- and fast-rise phases. In particular, 
we consider the transition between the two stages, using an analytic approximation to 
determine the onset of the fast-rise phase. 

4.1. Procedure 

Eruptions on the limb were chosen for their heightened contrast and limited projection 
effects compared to those on the disk. Bulk motions are inferred by tracing the leading edges 
in 304 A observations and fitting curves to the resultant height-time profiles. Figure 10 
illustrates this process for a single event. To begin, we select the linear slice that best 
characterizes the overall trajectory of the filament (Panel A). Additional slices offset by 
two degrees in either direction are identically processed for error estimation. Emission along 
these lines is binned to 300 pixels and interpolated onto a uniform grid, yielding a spatial 
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Start time = 2013/01/21 17:00:31 


Figure 10. An illustration of the kinematics procedure using event No. 0620. (a) Trajectories used for 
tracking. The black slice is selected manually, and the others are offset by 2° in either direction, (b) 
Height-time image for the black slice in Panel A. The dashed line indicates the time shown in Panel A, 
the solid curve is the fit result, and the asterisk denotes the fast-rise onset point, (c) Output of the Canny 
edge detection algorithm applied to Panel B. The red pixels are used as individual height measurements, 
(d) Fit to the red pixels in Panel C. The dotted lines indicate the onset of the fast-rise phase, and the 
colors correspond to results from the different slices in Panel A. (e) Velocity and (f) acceleration profiles 
for the height profile in Panel D. A movie showing the evolution of the eruption along the main slice is 
provided on the catalog website. 


resolution of ~1.5”, which varies slightly depending on the length of the slice, dictated 
by FOV position. The resultant light curves are stacked against subsequent observations 
to produce height-time images like the one shown in Panel B. To further improve signal- 
to-noise, averages of 10 frames were used for each column measurement. This averaging 
corresponds to an effective exposure time of 29 seconds over a period of 2 minutes. Shorter 
averages were used for fast eruptions where a 2-minute time resolution was insufficient. 

The height-time images are then further processed to improve contrast, which is typi¬ 
cally limited to multiplying each row by its height to boost signal far from the limb and 
thresholding the image above some multiple of its median value. Smoothing is applied for 
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Figure 11. Left: Three examples with varying goodness-of-fits. Right: Histogram of height uncertainties, 
which correspond to the values required for Xjj = 1.0 for the fits to Equation 1. 


particularly noisy cases, and base differencing (background subtraction) is used for events 
where the contrast was especially low. The Canny edge detection algorithm (Canny, 1986) 
is then employed to extract leading edges from the height-time images. Panel C shows the 
application of the Canny algorithm to the image in Panel B. Pixels highlighted in red, 
which form the uppermost edge, are used as the individual height measurements. These 
points are extracted automatically once the edge detection has been applied, but their 
time range must be selected manually. This consideration is particularly important for the 
start time because it effectively dehnes the onset of the slow-rise phase. 

4.2. Analytic Approximation 

Several curves have been used in previous studies to ht the low-corona height evolution 
of eruptive prominences and CMEs in general. When there are distinct phases, the initial 
slow-rise prohle is typically ht with a line, implying constant velocity (e.g. Sterling and 
Moore 2005), though at least one previous study has also employed a very slight constant- 
acceleration curve (Joshi and Srivastava, 2007). The fast-rise phase and rapid acceleration 
that initiates it has been approximated by exponential (e.g. Goff et al. 2005; Williams et al. 
2005), exponential-plus-initial-velocity (Alexander, Metcalf, and Nitta, 2002), constant- 
acceleration (e.g. Gilbert et al. 2000; Kundu et al. 2004), and cubic (Schrijver et a/., 2008) 
curves. Cheng et al. (2013) present an equation that combines a linear curve to treat the 
slow-rise phase and an exponential to treat the fast-rise after a time offset: 

h{t) = + Ci{t- to) + C2, (1) 

where h{t) is height, t is time, and r, to, cq, ci, and C 2 are free parameters. This model 
has the distinct advantage that both phases can be approximated by a single-function ht 
to the full dataset. It also provides a convenient method for determining the onset of the 
fast-rise phase, which can be dehned as the point at which the exponential component of 
the velocity equals the linear (i.e. the total velocity is twice the initial), given by: 


tonset = rln(ciT/co) + to 


( 2 ) 
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Table 3. Equation 1 h{t) fit variable statistics for 78 events, t is in 
seconds relative to the start of observations, and h is in arcsec along 
the trajectories shown in Figure 13. 


Variable 

Min 

Max 

Mean ± SD 

Median ± MAD 

T 

260 

19000 

2700 ± 2500 

2100 ± 1400 

to 

2600 

130000 

35000 ± 22000 

33000 ± 16000 

Co 

37 

1400 

240 ± 200 

180 ± 130 

Cl 

0.00028 

0.020 

0.0029 ± 0.0027 

0.0026 ± 0.0015 

C2 

30 

310 

120 ± 61 

110 ± 49 


We find this approximation to be best snited for onr stndy, which aims to characterize the 
the slow- and fast-rise phases of many events in a nniform way. Fitting is accomplished 
nsing MPFIT, a non-linear least sqnares cnrve htting package for IDL (Markwardt, 2009). 
Panels D, E, and F of Fignre 10 show the ht resnlt and its time derivatives for onr example 
event. Based on this approximation, we hnd that the initial slow-rise velocity is 1.8 ± 0.1 
km s"^, and the maxim um velocity in the AIA FOV is 91 ± 7 km s"^. The onset of the 
fast rise phase occurs at 00:48 UT ± 5 minutes at a distance of 55 ± 5 Mm along the ht 
trajectory, which begins with an initial height of 0.015 Rq (~10 Mm). The onset point is 
therefore 65 Mm above the limb along the fit trajectory, corresponding to a radial height 
of 64 ± 5 Mm. At the onset of the fast-rise phase, the acceleration is 0.57 ± 0.03 m s"^, 
and the hnal acceleration is 29 ± 4 m s"^. Table 3 lists the range of ht parameters found 
for the 78 events that can be satisfactorily described by this model. Details on the events 
that were not ht are given in §4.4. 

Two strategies are employed to quantify uncertainty, and errors quoted are the sum of 
both. The hrst is to identically process two adjacent slices ohset by 2° on either side of 
the original. On average, standard deviations from these results account for ~65% of the 
velocity, acceleration, and time uncertainties, but they may contribute as little as 30% or 
as much as 90% for a particular event. The second strategy is to perform 100 Monte Carlo 
(MC) simulations to estimate uncertainties from the ht parameters by randomly varying 
our height measurements within some assumed error and rehtting Equation 1. Since there 
are no standard errors for height measurements obtained by our edge detection method, 
the assumed height errors are chosen to yield a reduced chi-squared {xr) of 1-0 for the 
ht. This uncertainty is therefore a proxy for the goodness of ht. In the example above, 
height errors of 5.2" are required for Xr = 1-0; which is very close to the average value 
of 5.8" across all events. 5.8" corresponds to ~10 AIA pixels and around 4 pixels on the 
binned height-times images used to obtain the measurements. These values account for 
58% of the fast-rise onset height errors on average, while uncertainties associated with the 
three separate slices and MC realizations typically account for 28% and 14% of the onset 
height uncertainties, respectively. The left panel of Figure 11 shows the data, errors, and 
associated hts for three examples with varying height uncertainties, and the right panel 
shows the total distribution of height uncertainties. 

4.3. Solar Rotation 

As we will see in the next section, slow-rise velocities of < 1 km s"^ are found for a number 
of events and have also been reported in the literature (e.g. Sterling and Moore 2004; Isobe 
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Figure 12. The effect of solar rotation on apparent height (left) and velocity (right) over a 4-day period 
for a prominence with an actual height of 83 Mm, a latitude of 45°, and a Carrington longitude of -38.1°. 

et al. 2007; Liu et al 2012b; Regnier, Walsh, and Alexander 2011; Reeves, McCauley, and 
Tian 2015). It is important to note that solar rotation alone can produce apparent rise and 
fall velocities of this magnitude. Accounting for this effect requires an understanding of the 
projection geometry and the longitudinal extent of the hlament, consideration of which is 
beyond the scope of this work. The effect is minor and should be at least partially washed 
out in our statistics, but we nevertheless include this discussion to demonstrate that very 
small velocities for specific events should be treated with some caution. 

Foullon and Verwichte (2006) provide the following expression for the height of an 
arbitrarily thin prominence sheet above the limb, as seen from Earth: 

+ cos(L-Lo)-l, (3) 

where h{L) is apparent height, ho is actual height, L is the Carrington longitude of the 
limb, and Lq is the Carrington longitude of the prominence. Figure 12 illustrates this effect 
over a 4-day period for a hypothetical prominence with a height of 83 Mm (our average 
fast-rise onset height) at a latitude of 45° and a Carrington longitude (Lq) of -38.1°, which 
corresponds to that of the east limb at 2014/01/01 00:00 UT. We see that rise and fall 
velocities approaching 1 km s"^ can be observed when a prominence begins to rotate onto 
and off of the limb. This could cause non-negligible over- or underestimations of the slow- 
rise velocites for the slower events in our sample, depending on their projection geometries 
and longitudinal extents. 

4.4. Kinematics Results 

Figure 13 summarizes the kinematics results, showing the positions and velocity prohles 
of each eruption. Of the 106 events, 78 (74%) have distinct slow- and fast-rise phases 
that are well-ht by the linear-plus-exponential model proposed by Cheng et al. (2013). 
These include 44 QS, 22 PC, 8 IP, and 4 AR eruptions. Representative examples of 
the remaining events that cannot be satisfactorily described by Equation 1 are shown 
in Figure 14. Corresponding to the hgure panels, these events cannot be ht because they 
a) most commonly do not exhibit slow-rise phases or have no marked transition between 
phases, b) have rise profiles too complicated to fit into the simple two-phase dichotomy, c) 
have distinct phases but the transition into the fast-rise is decidedly not exponential, or 
d) exhibit dynamics that make it difficult to infer bulk motions from that of their leading 
edges. These events include 15 QS, 9 PC, 2 IP, and 2 AR events, and they are represented 
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Figure 13. Positions and velocity profiles from the kinematics study. The asterisks denote the fast-rise 
phase onset point for events fit with Equation 1. Events without onset points are described by exponential 
functions. The dashed circle represents the average radial onset height of 83 Mm, and the rounded corner 
edges mark the boundary of the telescope’s filter. 



23:06 00:13 01:20 02:26 17:13 20:00 22:46 01:33 

Start time = 2010/11/16 23:00:32 Start time = 2011/11/04 15:02:32 


Figure 14. Representative examples of the “bad” events that could not be described by Equation 1 and 
were excluded from further analyses. See §4.4 text for details. 
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Table 4. Kinematics Results. 


Parameter 

Min 

Max 

Mean ± SD 

Median ± MAD 

(/r±a)“ 

(AQ'’ 

^onset*^ (Mm) 

27 

237 

83 ± 41 

71 ± 32 

4.3 ± 0.51 

0.35 

^^onset*^ (Mm) 

2.6 

89 

29 ± 22 

22 ± 17 

3.1 ± 0.86 

0.37 

^^slow —rise (hrs) 

0.36 

23 

4.4 ± 3.7 

3.2 ± 2.7 

1.1 ± 0.93 

0.37 

'^slow-rise (km S ) 

0.21 

16 

2.1 ± 2.1 

1.9 ± 1.1 

0.50 ± 0.67 

0.81 

"^fov —max (km S ) 

6.0 

417 

106 ± 85 

82 ± 60 

4.4 ± 0.83 

0.18 

(km s'Q 

13 

753 

183 ± 169 

126 ± 110 

4.8 ± 0.84 

0.38 

flonset (m S'2) 

0.032 

29 

1.9 ± 3.7 

0.82 ± 1.5 

-0.28 ± 1.1 

0.53 

^fov —max (l^ ® ) 

0.29 

970 

111 ± 183 

40 ± 96 

3.6 ± 1.5 

0.38 

“500"'' (m s-2) 

0.64 

1760 

199 ± 362 

57 ± 175 

4.0 ± 1.6 

0.43 

f 

^onsef^ 

0.69 

2.0 

1.1 ± 0.28 

1.1 ± 0.22 

0.089 ± 0.25 

0.56 


“Lognormal fit parameters. 

^Anderson-Darling goodness-of-fit statistic; low values imply better fits. 
“Radial fast-rise onset height. 

“^Slow-rise phase displacement along fit trajectory (relative onset height). 
“Parameters at 500" above the limb; excludes 6 confined events. 

•^Decay index at the fast-rise onset point; see Section 5.1. 


in Figure 13 with exponential fits. While it would be interesting to consider these mishts 
in detail, they are excluded from all further analysis and left for future work. 

It is also clear that the kinematics study is biased toward QS and PC eruptions from 
the proportion of AR and IP events in this sample compared to the catalog. This selection 
effect is primarily due to the semi-automated method described in §4.1, which requires 
that the trajectory be free of additional structures that cannot be easily subtracted away. 
The more dynamic environments of AR and IP hlaments thus precluded a clean isolation 
of the eruption front by our edge detection method in a number of additional events 
that we attempted to process. This effect is exacerbated by the smaller size and lower 
initial heights of AR hlaments because they consequently spend less time above the limb. 
Several otherwise suitable events were also excluded because a signihcant portion of their 
slow-rise phases occurred behind the limb or on the disk. For this reason, our AR sample 
includes somewhat larger hlaments than average for that type. This bias combined with 
the small number of events means that our AR results may not accurately represent the 
general population, and additional work is therefore needed to obtain robust statistics. 
Rehnements to our technique can likely improve its applicability, but that is beyond the 
scope of this work. 

Figures 15, 16, and 17 show distributions of the fast-rise onset heights, slow-rise du¬ 
rations, velocities, and accelerations for the 78 events ht with Equation I. Cumulative 
distributions are also plotted for QS and PC events, along with a single curve for AR 
and IP eruptions, which are combined due to their small sample sizes. The histograms are 
best described by lognormal distributions, reminiscent of CME speeds (Yurchyshyn et al, 
2005; Zhang and Dere, 2006; Bein et a/., 2011), and lognormal hts are plotted in red. 
An independent variable, x, is said to be lognormally-distributed when ln(x) has a normal 
distribution. Note that Figure 17 plots the natural logarithm of acceleration because of the 
particularly large ranges and therefore exhibits a normal distribution. Physically, lognormal 
distributions may arise when a quantity is the multiplicative product of several independent 
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Table 5. Mean ± SD Kinematics Results by Type. 


Parameter 

All 

AR 

Type 

IP 

QS 

PC 

Event Count 

78 

4 

8 

44 

22 

^onset (Mm) 

83 ± 41 

48 ± 16 

60 ± 24 

82 ± 42 

101 ± 38 

^^onset (Mm) 

29 ± 22 

13 ± 11 

20 ± 9.4 

25 ± 19 

43 ± 24 

^^slow —rise (hrs) 

4.4 ± 3.7 

1.2 ± 0.82 

2.2 ± 1.9 

3.7 ± 2.8 

7.2 ± 4.4 

'^slow —rise (km S ) 

2.1 ± 2.1 

2.3 ± 0.41 

4.1 ± 4.8 

1.9 ± 1.5 

1.7 ± 0.89 

^fov —max (km S ) 

106 ± 85 

158 ± 95 

218 ± 111 

108 ± 74 

53 ± 33 

^^500" (km s'l) 

183 ± 169 

359 ± 180 

390 ± 216 

183 ± 146 

80 ± 52 

flonset (m S'2) 

1.9 ± 3.7 

4.1 ± 3.4 

6.6 ± 9.2 

1.4 ± 1.6 

0.60 ± 0.49 

*^fov —max (^ S ) 

111 ± 183 

320 ± 261 

336 ± 316 

97 ± 132 

21 ± 24 

« 500 " (m s-2) 

199 ± 362 

505 ± 402 

628 ± 589 

176 ± 305 

33 ± 47 

^onset 

1.1 ± 0.28 

1.0 ± 0.18 

1.1 ± 0.21 

1.2 ± 0.29 

1.1 ± 0.28 


variables, suggesting a multiplication of several independent processes, as opposed to the 
sum of independent variables for a normal distribution. 

The lognormal probability, /(x), and cumulative probability, F{x), distributions are 
given by: 


fix) = 


\/2'K( 


-exp 


TTUX 


(Inx — iJiY 

2^^ 


; Fix) = \ + ^erf 


( 


Inx 




V V2c 


(4) 


where (i is the mean and a is the standard deviation of ln(x). A least squares ht to the 
empirical cumulative probability distribution is used to approximate F{x), avoiding the 
unnecessary binning that would be needed to ht fix). We employ the modihed Anderson- 
Darling (A^) test to quantify the goodness-of-ht, which tells us if the empirical distribution 
is unlikely to be drawn from the htted distribution. The null hypothesis for this test 
is that the empirical distribution is drawn from the htted distribution, which can be 
rejected at the 0.05 signihcance level if > 0.752 and at the 0.25 level if > 0.470 
(D’Agostino and Stephens, 1986). Low values of A^ therefore imply better hts. Table 4 lists 
several statistical measures for each of the total distributions, including the lognormal ht 
parameters and A^ values, and Table 5 lists the arithmetic means and standard deviations 
for each quantity by hlament type. With reference to Bein et al. (2011) and Limpert, Stahel, 
and Abbt (2001), we use /r* = and a* = to refer to the median and multiplicative 
standard deviation of the htted distribution, respectively. The 68% conhdence interval for 
a lognormal distribution then spans from yi* -h a* to fi* x a*. 


4.4-1. Fast-Rise Phase Onset Heights & Slow-Rise Durations 

The left panel of Figure 15 shows the distribution of radial fast-rise onset heights (honset) 
derived from the prominence positions at tonset (Equation 2). We hnd mean and median 
values of 83 and 71 Mm, respectively, and the htted lognormal distribution has a median 
(/U*) of 74 Mm and a 68% conhdence interval that spans from 44 to 122 Mm. Compared to 
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Figure 15. Left: Radial fast-rise onset heights. Center: Slow-rise phase displacements along the fit tra¬ 
jectories. Right: Slow-rise phase durations. The histograms show the overall count distributions (left axes), 
the black curves show the cumulative probability distributions for specific types (right axes), and the red 
curves show lognormal fits to the histograms. 


QS eruptions, for which honset = 82 Mm, AR and IP hlaments tend to reach their critical 
points at 0.59 x and 0.74 x lower heights, respectively, while PC events transition into the 
fast-rise phase at 1.24x larger altitndes. We can also consider the distances traveled along 
the ht trajectories, which may be non-radial (see Fignre 13), by snbtracting the initial 
leading edge positions from their positions at fonset- This relative onset height (A/ionset) 
reflects the magnitnde of the slow-rise phase displacement, and its distribntion is shown in 
the middle panel of Figure 15. r* = 21.3 Mm for Ahonset, very close to the actual median 
value of 21.8 Mm, and the 68% conhdence interval ranges from 9 to 50 Mm. The relative 
onset heights are 0.51 x lower, 0.80x lower, and 1.75x higher for AR, IP, and PC events, 
respectively, compared to the QS average of 25 Mm. Some of the differences between both 
the radial and relative onset heights across hlament types can be attributed to a positive 
correlation between latitnde and onset height fonnd for the entire sample, which will be 
discnssed in §5.1 with respect to the decay index of the magnetic held. 

Liu et al. (2012a) also examine the fast-rise onset height, which they refer to as the 
critical height, nsing STEREO observations of 362 prominences that were detected and 
tracked by an automated system (Wang et al., 2010). Height-time measurements from the 
fast-rise phase were ht with a line, and the critical height was dehned as the intersection 
between this line and the height prohle prior to the fast-rise phase. This method would 
give somewhat larger heights if applied here, since our onset heights refer to the base of an 
exponential growth curve (Equation 2), bnt this diherence is very small compared to the 
range of heights fonnd. Lin et al. (2012a) hnd an average critical height of 77 Mm (0.11 
Rq), which is 6 Mm less than onr average but still well within one standard deviation. 

The right panel of Fignre 15 shows the distribntion of slow-rise phase dnrations, or the 
time reqnired to cover A/ionset- The average across all types is 4.4 honrs, and the type 
averages are sorted as expected, with AR and PC events at the low and high ends. The 
68% co nfi dence interval spans 1.2 to 8.0 honrs, which covers most of the values reported in 
the literature (e.g. Sterling and Moore 2005; Joshi and Srivastava 2007; Reeves, McCauley, 
and Tian 2015). The minimum slow-rise phases were 22 and 24-minntes for an IP and AR 
event, respectively. The short end fonnd in the literatnre is aronnd 10 minntes, reported by 
Sterling and Moore (2005) for an AR ernption, thongh the anthors note that the filament 
appeared to begin rising at least 10 additional minntes earlier. Chifor et al. (2006) report 
and AR hlament slow-rise dnration of 34 minntes for a separate event. Onr sample of AR 
events is qnite small for reasons described above, which likely contribntes to the paucity of 
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Initial Velocity (km s"') Peak Velocity in FOV (km s"') Velocity at 500” (km s"') 

Figure 16. Left: Slow-rise velocity. Center: Peak velocity in FOV. Right: Velocity at 500” above the limb 
(projected). The histograms, curves, and axes are as in Fig 15. 


brief slow-rise phases and biases our lognormal fit toward longer durations. AR eruptions 
may also exhibit only exponential growth (i.e. have no slow-rise phase), such as the two 
AR events that could not be ht by Equation 1 in our initial sample and the event described 
by Williams et al. (2005). In these cases, the slow-rise phase may be non-existent or may 
simply be indiscernible by the analysis technique or instrument resolution. At the other 
end of the spectrum, the longest slow-rise phase is 23 hours for a PC eruption. This event is 
an extreme outlier at 3 hours outside of the 95% conhdence interval upper bound (p* x a*‘^) 
for the lognormal ht, already expected to be biased toward long durations, and indeed we 
could hnd no counterparts in the literature. Liu et al. (2012b) report a slow-rise phase 
duration of ~10 hours, which is one of the largest values previously reported and is much 
closer to the longer slow-rise phases we observe, excluding the outlier. 

4 . 4 . 2 . Velocity & Acceleration Distributions 

Moving to the velocity distributions plotted in Figure 16, we hnd that the initial slow-rise 
speeds range from a minimum of 0.2 to a maximum of 16 km s"^. Recalling §4.2, these values 
are derived from the linear component of Equation 1 and dominate the total velocity until 
the exponential component of the curve contributes equally at fonset (Equation 2). The 
lognormal ht to the Ugiow-rise distribution is particularly bad, with > 0.752, meaning 
that we can reject the hypothesis that the velocities are lognormally-distributed with a 
p-value < 0.05. Instead, the histogram is double-peaked with enhanced slow-rise velocities 
around 1 and 2 km s"^. However, if we consider only QS events, the values follow a lognormal 
distribution very closely to = 0.25. (In fact, all of the QS-only distributions except 
for ffov-max somewhat better-ht by lognormal prohles than the corresponding total 
distributions, with A^ = 0.31 vs. 0.42. While we do not show the type-specihc histograms 
and ht prohles, the characters of their distributions may be gleaned from the cumulative 
distributions in Figures 15-17.) 

The total Ugiow-rise distribution deviates from log-normality primarily because the PC- 
only distribution, which has a large enough sample to draw meaningful conclusions from, 
is not lognormal to a signihcance of over 99% (A^ = 1.1). Instead, it is better ht by a single 
Gaussian with a mean of 1.7 and a variance of 1 km s"^ (A^ = 0.55). The low number of 
events in the AR and IP distributions preclude ns from concluding anything about their 
general fsiow-rise distributions, but the limited information hints at log-normality for IP 
and normality for AR eruptions. In light of this, the other panels of Figure 16 reveal an 
intriguing development. The middle panel shows the maximum velocities attained in the 
AIA FOV, and the right panel shows the ht velocity at a hxed height of 500" above the 
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Onset Accel. [In(m s'^)] Peak Accel. In FOV [ln(m s'^)] Accel, at 500” [ln(m s'^)] 
Figure 17. Left: Fast-rise onset acceleration. Center: Peak acceleration in FOV. Right: Acceleration 
500” above the limb (projected). The histograms, curves, and axes are as in Fig 15. 


at 


limb, which is roughly the maximum height that AIA can observe for structures near the 
FOV corners. These latter velocities are included because the maximum observable height, 
where ffov-max is generally taken, varies depending on FOV position (see Figure 13). Low- 
latitude hts are therefore projected for the 500" comparison, and 6 conhned events are 
excluded because the hlament material does not reach that height. 

Both of these velocity distributions exhibit robust lognormal prohles despite having 
emerged from a slow-rise distribution that is signihcantly less so. This is true even for 
the PC-only distributions (A^ 0.33). Lognormal prohles are also found for each of the 

acceleration distributions shown in Figure 17, which is particularly signihcant for the fast- 
rise onset accelerations (uonset) because it indicates that the velocity distribution shift 
begins immediately at the start of fast-rise phase. We might draw two inferences from 
this. First, the normal distributions of PC and perhaps AR slow-rise velocities hint at a 
more uniform process than the IP and QS counterparts. This hnding may be understood 
in terms of the more varied environments and size-scales of IP and QS hlaments, which 
provide a more diverse set of inhuences to impact the slow-rise pace. Second, the lognormal 
distributions for the fast-rise velocities of all types hint at the multiplicative contribution 
of several processes. This result may be understood in terms of the various, and not 
necessarily mutually exclusive, fast-rise onset mechanisms (to be discussed further in §5) in 
combination with environmental influences. We also see that the lognormal distribution of 
CME speeds reported in previous works emerges early in the height evolution of prominence 
eruptions. 

Now we turn briefly to how our velocity and acceleration values differ by hlament type. 
As with the CME speeds in §3.8, we hnd that IP eruptions are generally the fastest of the 
four types. The average IP slow-rise velocity is 4.1 km s"^, which is about I.8x, 2.2x, and 
2.4X faster than the AR, QS, and PC averages, respectively. Similar relationships between 
types are found for the other velocity and acceleration quantities. Accelerations at fonset 
range from a minimum of 0.03 m s"^ for a PC eruption to a maximum of 29 m s"^ for an 
IP event, with an overall average value of 1.9 and a 68% conhdence interval of 0.25 to 2.3 
m s"^. After the fast-rise onset, the velocities grow exponentially to reach values at 500" 
above the limb of between 13 km s"^ for the slowest PC eruption and 750 km s"^ for the 
fastest IP event. The median value upon exiting the AIA FOV is 82 km s"^. See Table 4 
and Table 5 for the full list of statistical parameters. These values are broadly consistent 
with those reported in the literature, much of which has already been referenced. 
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Figure 18. Latitude versus initial height (left), onset height (honset, middle), and relative onset height 
(Ahonset, right), where the latitudes are those of the onset points in Figure 13. 

4.5. Online Catalog Kinematics Content 

The online catalog inclndes a version of Fignre 10 and the associated movie for each of the 
106 events in the kinematics stndy. Text hies hies containing the height vs. time data, ht 
parameters, and essential resnlts are also provided. Information on hnding and nsing these 
hies is provided on the catalog website. 

5. Discussion 

A few interesting patterns emerge from our results that are relevant to the broader dis¬ 
cussion on hlament eruption triggering. Initiation mechanisms for the fast-rise, rapid 
acceleration phase (or for the entire eruption if there is no slow-rise) are generally divided 
between ideal MHD instabilities and fast magnetic reconnection processes (see review by 
Fan 2015). The former includes the helical kink instability, discussed in §3.5, and the 
torus instability, which occurs when a toroidal hux rope rises to a critical height above 
which it cannot be contained by the overlying magnetic held (Kliem and Torok, 2006). 
Considering reconnection triggers, hux emergence (Chen and Shibata, 2000) or cancellation 
(van Ballegooijen and Martens, 1989) may lead to reconnection that stimulates the erup¬ 
tion, or breakout reconnection in the overlying held may prompt an eruption (Antiochos, 
DeVore, and Klimchuk, 1999). These mechanisms are not necessarily mutually exclusive, 
and reconnection may also facilitate an eruption triggered by the catastrophic loss of MHD 
equilibrum (Lin and Forbes, 2000). Several mechanisms have also been proposed for the 
slow-rise phase, perhaps the most common of which is hux cancellation that increases the 
twist of a hlament, slowly buoying it until an MHD instability or reconnection process 
triggers rapid acceleration (e.g. Amari et al. 2010; Aulanier et al. 2010). Alternatively, 
Reeves et al. (2010) simulate the initial slow rise of a hux rope CME by means of resistive 
dihusion, and Fan and Gibson (2007) discuss the possibility of a hux rope emerging from 
beneath the photosphere. 

5.1. Onset Height, Latitude, &; Decay Index 

Figure 18 plots the initial prominence height, onset height, and relative onset height versus 
latitude, where the latitudes are those of the onset points in Figure 13. Recalling §4.4.1, 
the onset height (honset) is a radial height above the limb, while the relative onset height 
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Onset Height Decay Index 


Figure 19. Left: Average decay index between 42 and 105 Mm vs. latitude. Symbols and colors are as in 
Fig. 18. Right: Distributions of the decay indexes at the fast-rise onset height. The histogram, curves, and 
axes are as in Fig. 15—17. 


(A/ionset) refers to a displacement along a particular fit trajectory. We find a positive 
correlation between latitude and onset height, with a Pearson’s correlation coefficient (r) 
of 0.30 and a p-value of 0.008. The relationship strengthens if we instead consider the 
relative onset height (r = 0.43, p = 7x10“^) or slow-rise phase duration (r = 0.50, p = 
3x10“®). This effect is distinctly global in that if we divide the data into latitude bins of 
20°, for instance, the correlation within a given bin is insignihcant or non-existent. It is 
also distinct from a relationship between latitude and the initial prominence height, for 
which there is a weak and statistically insignihcant positive correlation (r = 0.14, p = 0.13) 
that disappears if IP and AR events are excluded. 

The maximum height of non-eruptive prominences was suggested as a parameter of 
the magnetic held by Leroy, Bommier, and Sahal-Brechot (1984). Makarov et al. (1992) 
expanded on this, using height variations among prominences of the same latitude as 
a proxy for longitudinal variations in the background held strength. Makarov (1994) also 
reported that the height of polar crown prominences (latitudes greater than 55°) decreased 
with increasing latitude using observations of many prominences over several solar cycles. 
We see this ehect to some extent in our initial height data for the PC-only sample, which 
exhibit a minor negative correlation (r = -0.21, p = 0.35) between latitude and initial 
height. Upon erupting, the PC-only data shift to no correlation with latitude for onset 
height and a minor positive correlation (r = 0.17, p = 0.45) for relative onset height. 

The critical height of a flux rope, above which no stable equilibria exist, can be related to 
a critical threshold in the vertical gradient of the magnetic held, or decay index (Filippov 
and Den, 2000): 


d(lnB) 
(i(ln/i) ’ 


( 5 ) 


where B is the strapping held and h is height above the photosphere. This is the essential 
concept of the previously-mentioned torus instability, which refers specihcally to toroidal 
(curved) hux ropes but may sometimes be used more generally to encompass all geometries. 
The critical decay index (ucr) for loss of conhnement depends principally on the assumed 
hux rope conhguration and is canonically 1.0 for the straight (van Tend and Kuperus, 
1978) and 1.5 for the toroidal (Bateman, 1978) thin current channel approximations. Both 
geometries represent limiting cases of the same physical process (Demoulin and Aulanier, 
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2010). Recent theoretical calculations and numerical simulations have found Ucr values 
between 1.0 and 2.0 under various assumptions that better approximate coronal conditions 
(Kliem and Torok, 2006; Fan and Gibson, 2007; Aulanier et ai, 2010; Demoulin and 
Aulanier, 2010). If we assume that many of the events in our sample are driven to erupt 
by this sort of mechanism, then a positive correlation between latitude and fast-rise onset 
height implies that the decay index is smaller at high latitudes for a given altitude. 

To test this hypothesis, we estimate the decay index along our eruption trajectories 
using potential held source surface (PFSS) extrapolations of the coronal magnetic held 
based on line-of-sight magnetograms from the SDO’s Helioseismic and Magnetic Imager 
(HMI, Scherrer et al. 2012). Decay index maps are obtained for the plane of the sky through 
Sun-center using the SolarSoft IDL software package FORWARD® (Gibson et al, 2014), 
which is capable of forward modeling various observables but in this case simply provides 
a convenient interface for SolarSoft’s PFSS package®. Details on the PFSS model are given 
by Schrijver and De Rosa (2003). For each event, we record the decay index at the fast- 
rise onset height (uonset) and the average decay index between 42 and 105 Mm above the 
limb (uavg), the range of which is chosen for consistency with Liu (2008). The left panel 
of Figure 19 shows Uavg versus latitude, and we do hnd the expected negative correlation 
(r = -0.22, p = 0.06). This result, combined with the slight positive correlation between 
initial height and latitude, nicely explains the positive correlations between latitude and 
the related measures of fast-rise onset height, relative onset height, and slow-rise duration. 

The right panel of Figure 19 shows a histogram of the fast-rise onset height decay indexes, 
and basic statistics are listed in Tables 4 and 5. It is interesting to note that the Uonset 
distribution clusters around 1.0, closer to the straight flux rope approximation than the 
curved. The average value is 1.1 and is essentially the same for all four hlament types, 
despite their differences in kinematics parameters, particularly onset height. This result 
is generally consistent with Filippov and Den (2001) and Filippov (2013), who also hnd 
critical decay indexes of around 1.0 for observed quiescent hlament eruptions and therefore 
suggest that the forces associated with the axis of curvature are essentially unimportant for 
the equilibria of coronal hux ropes. However, we also hnd that a signihcant fraction (27%) 
of our events have critical decay indexes greater than 1.3. Indeed, our results span the full 
range of aforementioned theoretical Ucr values derived from various hux rope and simulation 
conhgurations. We also hnd that 38% of our events begin their fast-rise phases with decay 
indexes of less than 1.0 and 15% do so with Uonset < 0.9. This hnding suggests that an 
appreciable fraction of the eruptions are initiated by other means, such as reconnection or 
the kink instability. 

Some limitations to our decay index estimates should be noted. First, we are looking 
at regions on the limb, where the PFSS extrapolations are constrained by observations 
a few days before or after the eruptions for the west and east limbs, respectively. Such 
time ohsets are especially limiting for active regions, which are very likely to have evolved 
over those periods. Second, we use single values in the plane of the sky, which may not be 
perfectly aligned with a given hlament. These considerations are particularly important for 
the discussion of eruption mechanisms because relatively small changes in the decay index 
can imply signihcant diherences in the initiation process. They might also explain two 
peculiarities in our decay index results. The curvature of AR region hlaments, with their 
shorter lengths, is generally thought to be more important compared to that of their more 


®SolarSoft IDL FORWARD Package: http://www.hao.ucar.edu/FORWARD/ 
^SolarSoft IDL PFSS Package: http://www.lmsatcom/~derosa/pfsspack/ 
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Figure 20. Empirical cumulative probability distributions for events with and without apparent twisting 
motions (see §3.4). The right panel (CME speed) covers the entire catalog, while the others include the 78 
events in the kinematics study fit with Equation 1. 


extended QS counterparts, meaning that AR Ud- thresholds should be closer to the toroidal 
approximation of ~1.5. We instead hnd an average Uonset of 1-0 for the 4 AR events in our 
sample, which might be due to the PFSS limitations. There is also a conspicuous outlier 
in Figure 19, with a high latitude and large decay index. This hlament (No. 0934) appears 
to be inclined such that the plane-of-sky decay index is not appropriate. More careful 
modeling and the use of STEREO observations to triangulate the hlament positions would 
improve the decay index estimates. We will discuss possibilities for additional work on this 
topic in §5.3. 

5.2. Kinematics &; Twist 

It is also interesting to consider possible relationships between kinematics and apparent 
twisting motions. Referring to §3.4, “apparent twist” may include motion that is induced 
during the eruption by surrounding features (i.e. the “roll effect”) or, more commonly, 
motion that is likely to be indicative of the intrinsic hux rope twist, such as untwisting 
motions observed in hlament footpoints near the end of an event. We hypothesize that 
on average, eruptions exhibiting twisting motions will have had greater helicity at their 
fast-rise onset points compared to those with no apparent twisting motions. Given this 
expectation, it is interesting to note that all but 2 of the 18 events with no apparent twist 
lie very close to or above the ht lines in the middle and right of Figure 18 (latitude vs. 
onset-height). On average, events with no twist begin their fast-rise phases at 1.4x larger 
heights, and comparing the cumulative distributions in the left panel of Figure 20 indicates 
a statistically signihcant difference (A| = 3.0 and p = 0.02). This result may be attributed 
to differences in the eruption mechanisms that preferentially affect more or less twisted 
hlaments. One possibility is that highly-twisted flux ropes are more susceptible to the kink 
instability compared to less-twisted structures, which are more likely to succumb to the 
torus instability. These mechanisms should produce populations of lower and higher fast- 
rise onset heights, respectively, which would be superposed on a scatter of reconnection- 
driven events. The right panel of Figure 18, in particular, exhibits a pattern somewhat 
similar to this. 

Fan and Gibson (2007) hnd from MHD simulations that the held decreases with height 
more slowly for the kink instability compared to the torus instability. We might therefore 
expect to see diherences in the distributions of fast-rise onset height decay indexes (uonset) 
between our two twist groups if there are indeed systematic diherences in their initiation 
mechanisms. The second panel of Figure 20 shows that there is no such diherence, and the 
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same is true if we instead consider the average decay index between 42 and 105 Mm. This 
may be because of the limitations to our decay index estimates described in the previous 
section or because our sample is too small to detect a slight difference. Interestingly, Liu 
(2008) hnds there to be no signihcant difference between the average decay indexes of 
kink and torus instability events, in contrast to the MHD simulations, and limitations 
to the PFSS modeling are also cited as a possible reason for this discrepancy. The fact 
that we hnd no difference in the decay indexes between events with and without apparent 
twisting motions therefore does not preclude differences in their initiation mechanisms. 
Moreover, we do not expect that all events in the twisted sample are initiated by the kink 
instability-most probably are not. Rather, we suggest simply that the likelihood is greater 
than for events without apparent twisting motions, which might explain the difference in 
onset height distributions. It may also be that more highly twisted hlaments are likelier 
to facilitate a reconnection process at some point during their slow-rise phases before the 
torus instability threshold is reached. 

The velocities plotted in the three rightmost panels of Figure 20 also tell an interesting 
story. We see from the middle panel that events with and without apparent twisting motions 
have no general differences in their slow-rise velocities. If, as we hypothesize, the twisting 
motions indicate greater helicity on average, then this similarity suggests that the flux rope 
helicity upon eruption is largely decoupled from the slow-rise velocity. This constraint is 
useful for modeling efforts, and it might lead one to speculate that there is a comparative 
homogeneity in the slow-rise mechanisms across events with disparate fast-rise mechanisms. 
As with the onset heights, velocity and acceleration differences emerge between the twist 
groups during the fast rise phase. We hnd that events with apparent twisting motions 
compared to those without have hnal speeds in our kinematics study that are 1.3 x faster 
and CME speeds across the entire catalog that are 1.2 x faster on average. While these 
enhancements are not very large, the cumulative distributions in Figure 20 indicate that 
the differences are statistically signihcant, with A| = 1.9 & p = 0.05 for the AIA speeds 
and A| = 3.4 Sz p = O.OI for the CME speeds. This suggests that hux rope helicity may be 
positively correlated with CME speed, though a quantitative measure of helicity is needed 
to make this claim. 

Two previous studies are particularly relevant to this possibility. Sung et al. (2009) found 
a positive correlation between the square of the initial CME speed and the magnetic helicity 
per unit length of the corresponding magnetic clouds at I AU for 34 events. Park et al. 
(2012) examined the average helicity injection (H) of CME-producing active regions for a 
few days leading up to 47 eruptions using Michelson Doppler Imager (MDI) magnetograms 
and found a positive correlation between H and CME speed. This result is consistent with 
our hnding, though it relates specihcally to active region CMEs, while we hnd roughly the 
same enhancement across all hlament types. However, when divided by type, our result 
remains statistically signihcant only for QS and PC events (A| = 3.2 & p = 0.02), for 
which there are over twice as many events compared the AR and IP sample (A^ = -0.12 
& p = 0.39). We also note from §3.4 that conhned eruptions are signihcantly less likely 
to exhibit twisting motions (44% vs. 60% for the general population), implying lesser 
helicity on average. Conhned events are often arrested by the overlying held because there 
is insufficient energy released to push through, and the fact that we hnd signihcantly lower 
rates of apparent twisting motions in these eruptions suggests, albeit somewhat tenuously, 
that we are reasonable in considering the twisted sample to be generally more energetic. 
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5.3. Future Work 

A recently-submitted paper by Su et al. (2015) presents detailed modeling and observations 
for one of the polar crown events in our sample, which includes a more careful estimate 
of the decay index that employs STEREO observations to better constrain the location of 
the hlament. This study implicates the torus instability as the driving mechanism through 
a comparison of the decay index to our fast-rise onset height. Reeves, McCauley, and Tian 
(2015) present a related study on a prominence eruption observed by the Interface Region 
Imaging Spectrograph (IRIS) that includes a comparison of the eruption kinematics to 
brightenings observed by both IRIS and AIA, which suggest that reconnection below the 
prominence triggered the fast-rise phase. The techniques employed by these works can be 
readily applied to our kinematics sample to estimate relative frequencies of the fast-rise 
onset mechanisms. The suggestive relationship between apparent twisting motions and 
CME speed can also be explored by quantifying the twist or overall helicity in quiescent 
hlaments. This might be done using helicity injection measurements similar to that of 
Park et al. (2012) applied to HMI observations of events in the catalog observed on the 
disk or west limb, though the low flux densities in quiet sun regions may preclude reliable 
measurements. Alternatively, observations of untwisting motions similar to those presented 
by Yan et al. (2014) might be used to measure the extent to which hlaments are twisted 
and then to compare with CME speed measurements. 

6. Conclusion 

This paper details the hrst statistical study of prominence eruptions observed by the 
Atmospheric Imaging Assembly (AIA) aboard the Solar Dynamics Observatory (SDO). 
An online catalog containing 904 eruptions was developed, listing several properties for 
each event that are intended to aid the community in identifying eruptions for future 
work (§2). Filaments are typed based on environment (active region (AR), intermediate 
(IP), quiescent (QS), or polar crown (PC)) along with the extent to which they erupt 
(full, partial, or conhned), and comparisons between the type groups are made throughout 
(§3.1). A sample of 106 limb events was drawn from the catalog for a kinematics study 
that characterizes the distinct slow- and fast-rise phases often exhibited by prominence 
eruptions (§4). Our main results are summarized below: 

1. Symmetric eruptions are found to be somewhat more common than asymmetric (48% 
vs. 38%), and symmetric events are 1.5 x more likely to be full eruptions due to enhanced 
mass draining and detachments in asymmetric events (§3.2). Relatively few eruptions 
exhibit non-radial (12%) and sideways (11%) trajectories. Sideways trajectories are more 
commonly found in AR and IP events, which are likelier to be deflected strongly away 
from ARs. Sideways events are also more frequently conhned or partial eruptions, likely 
because of interactions with neighboring structures that may arrest the eruption or 
facilitate additional draining (§3.3). 

2. Vertical threads, perpendicular to the hlament spines, are observed for 58% of events, 
which is interesting given their perplexing nature. The fraction of events with vertical 
threads by type scales with the characteristic sizes of the four types, with 11% for AR and 
74% for PC eruptions. We note that vertical threads appear in ARs only for particularly 
large hlaments, often in decaying active regions with well-separate polarities or at the 
boundary of what might instead be considered an intermediate hlament (§3.6). Coronal 
cavities are also found most frequently in PC events, as expected. Overall, cavities are 
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absent in 49% of limb events, present in 34%, and become visible during the eruptions 
in 17% (§3.7). 

3. 72% of our filament eruptions are associated with white-light CMEs observed by LASCO, 
and similar rates are found for each of the hlament types. The CME speeds associated 
with AR and IP events are generally faster than those associated with QS and PC events 
by a factor of 1.3. IP-CMEs tend to be 1.1 x faster than AR-CMEs, but this difference 
is not statistically signihcant (§3.8). 

4. 106 limb events were selected from the catalog for a kinematics study, including 59 QS, 
31 PC, 11 IP, and 6 AR prominence eruptions. Height-time plots were constructed for 
each event using 304 A observations, and the Canny edge detection algorithm was used 
to extract individual height measurements (§4.1). A majority (74%) of the rise prohles 
are well-ht by an approximation introduced by Cheng et al. (2013), which combines a 
linear equation to treat the slow-rise phase and an exponential to treat the fast-rise. The 
onset of the fast-rise phase is dehned as the point at which the exponential component of 
the curve equals the linear (§4.2). We also illustrate that apparent rise and fall velocities 
approaching 1 km s"^ can be observed due to solar rotation alone (§4.3). 

5. The average values for the radial onset height (honset), relative onset height (Ahonset), 
and slow-rise (Atgiow-rise) phase duration are 83 Mm, 29 Mm, and 4.4 hours, respec¬ 
tively. As expected, the lowest heights and shortest slow-rise phases are found for AR 
eruptions, and the largest heights and longest slow-rise phases are found for PC erup¬ 
tions (§4.4.1). The average values for the slow-rise velocity (fsiow-rise), maximum FOV 
velocity (ffov-max), and velocity at 500" above the limb (usoo") are 2.1, 106, and 183 
km s"^, respectively. The average onset acceleration (uonset), maximum FOV acceleration 
(ofov-max), and acceleration at 500" (asoo") are 1.9, 111, and 199 m s"^, respectively. 
As with CME speed, we hnd that IP events are generally the fastest of the four types 
(§4.4.2). See Table 4 for the full list of statistical parameters and Table 5 for type 
comparisons. 

6. The kinematics parameter distributions are best described by lognormal probability 
distributions similar to that of CME speeds, indicating that this pattern emerges in 
the low-corona (§4.4). An exception to this is the slow-rise phase velocity distribution 
for PC and possibly AR eruptions, which follow Gaussian distributions that shift to 
log-normality during the fast-rise phase. This may suggest a uniformity in the slow- 
rise process for PC and possibly AR hlament eruptions compared to their QS and IP 
counterparts, whose more varied environments and size scales may provide more diverse 
inhuences (§4.4.2). 

7. We hnd a positive correlation between latitude and fast-rise onset height, with correla¬ 
tion coefficients (r) of 0.30 and 0.43 for the radial and relative onset heights, respectively. 
A corresponding correlation between latitude and slow-rise phase duration is also found 
(r = 0.50). This is a global ehect that diminishes if the eruptions are binned into latitude 
groups and is therefore likely to be a product of the global magnetic held. We interpret 
these correlations in terms of the decay index of the vertical magnetic held strength, 
which we hnd to exhibit a negative correlation with latitude using average decay indexes 
obtained from PFSS extrapolations between heights of 42 and 105 Mm. High-latitude 
events thus tend to have larger onset heights and longer slow-rise phases because they 
tend to have smaller decay indexes at a given height. We also hnd that the distribution 
of decay indexes at the onset of the fast-rise phase spans the full range of theoretical 
torus instability or loss of conhnement critical thresholds (~l-2), with an average value 
(1.1) consistent with the straight current channel approximation. A number of events 


SOLA: filament_catalog_paper.tex; 11 May 2015; 0:17; p. 33 



McCauley et al. 


also exhibit critical decay indexes of < 1.0, indicating initiation by other means, snch 
as reconnection or the kink instability (§5.1). 

8. A majority of events (60%) exhibit apparent twisting motions, and IP hlaments have the 
highest rate (68%). This inclndes several types of twist, of which nntwisting motions seen 
in hlament legs in the latter stages of an ernption are most common. We hypothesize 
that events with twisting motions have greater helicity on average, which is important to 
onr interpretation of the relationships found between twist and eruption kinematics (see 
below). Conhned eruptions exhibit twisting motions at a much lower rate (44%), which 
we interpret as a consequence of their comparatively limited magnetic energy (§3.4). A 
small fraction of events exhibit writhing motions suggestive of the kink instability (16%), 
but half of these are labeled as ambiguous. As expected, writhed events are much more 
likely to be partial or conhned than the general population (§3.5). 

9. Statistically signihcant differences are found between events that exhibit signs of twist 
and those that do not. Twisted events transition into their fast-rise phases at lower 
heights and have faster low-corona speeds, which might rehect different populations 
of fast-rise onset mechanisms that preferentially affect more or less twisted hux ropes. 
We also hnd that events with apparent twisting motions have faster coronagraph CME 
speeds by a factor of 1.2 across all events in the catalog. Though not very large, this 
difference is statistically signihcant {p = 0.01). Park et al. (2012) showed that helicity 
injection is positively correlated with CME speed for AR events, and our results indicate 
the same relationship for hlament eruptions of all types. We also hnd that there is no 
diherence in slow-rise speed for events with and without apparent twisting motions, 
which may imply that the hux rope helicity upon eruption is largely decoupled from the 
slow-rise mechanism and that there may be a comparative homogeneity in the slow-rise 
mechanisms across events with disparate fast-rise mechanisms (§5.2). 
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